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ABSTRACT 

The effect of slip condition on unsteady flow of an electrically conducting fluid trough a channel filled with 
saturated porous medium in the presence of transverse magnetic field and radiative heat transfer and mass transfer is 
studied. The dimensionless governing equations are solved using oscillatory flow conditions. The results are obtained for 
velocity, Sherwood Number and shearing stress for different parameters like Schmidt number, time, magnetic parameter, 
Darcy number, Renold number, Peclet number, rarefaction parameter etc. The flow characteristics are discussed and shown 
by means of graphs and tables. 
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INTRODUCTION 

It is well known that at the microscopic level the fluid velocity matches the velocity of solid boundary. While in 
no slip boundary condition has been proven experimentally to be accurate for a number of microscopic flows, it is an 
assumption but based on physical principles. About two hundreds years ago Navier proposed a general boundary 
conditions for the fluid slip at solid boundary. That velocity of solid surface is proportional to shear stress at the surface[l] 
.i.e 

,dV x 
V = h — - 
* dy 

where h is the slip co-efficient, if ft = Q then no slip condition obtain, if k (finite)fluid slip occurs at the 
wall. The above relation is linear but one can establish a nonlinear relationship of the slip flow. 

The fluid slippage phenomena of the solid boundaries appear in many applications such as micro channels or 
nanochannels and in application when surface is coated with special coatings such as thick monolayer of hydrophobic 
octadecyltrichlorosilance [9]. Recently , several researchers have suggested that the no slip boundary condition may not 
suitable for hydrophilic flows over hydrophobic boundaries at both micro and nano scale[4,7-8].The effect of fluid slippage 
at the wall for coutte flow are considered by Marques et al.[6]under steady state condition and only for gases. The effect of 
slip condition on MHD steady flow in a channel with permeable boundaries has been discussed by Makinde and 
Osalusi[ll]. Oscillatory free convection flow of a second order fluid from vertical plate have been studied by Mishra and 
Panda [3]. Soundalger et al. [5] have studied the effect of mass transfer in the flow past an infinite vertical oscillating plate 
in presence of constant heat flux .Makinde and Mhone[10] have studied the heat transfer on MHD oscillatory flow in a 
channel filled with porous medium. Senapati et.al.[12] have studied the effect of the mass and heat transfer effect on MHD 
oscillatory flow in a channel filled with porous medium. 
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In the recent paper we investigate the Effect of slip condition on unsteady MHD oscillatory flow in a channel 
filled with porous medium with heat radiation and mass transfer which is an extension work of Senapati et.al [12] 

FORMULATION OF THE PROBLEM 

The flow of a conducting optically thin fluid in a channel of width 'L' filled with saturated porous medium under 
the influence of an external applied homogenous magnetic field and radiative heat transfer is considered. It is assumed that 
Dufour and Soret effect viscous and joulean dissipation have been neglected. X -axis is taken along one wall of the channel 
and Y -axis is normal to it. 

A uniform magnetic field of strength H has been applied perpendicular to the wall of the channel and the 
magnetic permeability |i e is constant throughout the field. T and T w are the temperature and C and C w are concentration 
at both walls of the channel. Then by usual Boussinesq's approximation the unsteady flow is governed by following 
equations: 

dt p ox Qy p K 

„ dt d 2 T da 

pC — - = k — (2) 

P dt ay 2 dy 

dC d 2 C 

= D —^ (3) 



dt dy 



Following Cogley et al. [2], it is assumed that the above plate is heated at a constant temperature so the radiative 
heat flux is given by 

^ = 4« 2 (r -r) ( 4) 

dy 

where u is the axial velocity, P is the pressure, g is the gravitational force is radioactive heat flux, (3 is the co- 
efficient of expansion due to temperature, (3 C is the co-efficient of expansion due to concentration, C is the specific 

heat at constant pressure, k is thermal conductivity, K is permeability co-efficient of porous medium,B (|i e H o ) is the 
electromagnetic induction, <J is the conductivity of fluid, p is fluid density and V is the kinematic viscosity co-efficient 

The boundary conditions are as follows: 



, du . 
" =h —,T =T W ,C =C w fory =L 
dy 



u =0,T =T ,C =C Q fory =0 

Let us introduce the dimensionless quantities as: 



(5) 
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, = £^., G , = 'ftfc- c ^,P, = 



C, - C 



u 



pv 



2 T 2 



,N 2 = 



4a L 



— ,? e =R e .? r ,D a =^,h = - 
k L 2 L 



(6) 



J 



where D is mass diffusion, G r is Grashof number , G m modified Grashof number, H is Hartmann number , P r is 

Prandd number, S c is Schmidt number, P e is Peclet number , D a is Darcy number , N is radiation parameter, U is flow 
mean velocity, h is rarefaction parameter and R e is Renold number. 

With the help of boundary conditions (5) and equation (4) the equations (1) to (3) reduce to 

du <9P d u 



e dt dx by 2 m v r 

dt BY 2 
dl = d^_ 

e c dt by 2 
1 1 

Where S = 

D a 

The boundary conditions in non -dimensional form are 

u = h — ,0 = 1, <j> = lfory = 1 

dy 

u=O,0 = O,0 = Ofory = O 

METHOD OF SOLUTION 

As the flow is purely oscillatory, in order to solve the equations (7) to (9) Jet us introduce 



dx 



= Ae iM ,u(y,t) = u o ( y y M ,0(y,t) = o (y)e iM ,<p(y,t) = <P o (y)e i 



where 0) is the frequency of oscillation and A. is constant. 
Then equations (7) to (9) becomes 

d ^ U f ~ m i 2 "o + G r Q + GJ Q +1 = 
dy 

d 2 n 2 „ „ 
— f + m 2 O =0 



dy 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



(13) 
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dy 2 



-m 3 =0 



where m, = ^S 2 +H 2 + iR e O) , m 2 = ^/n 2 -P t 
By solving the equations (12) to (14) we get 
sinh(m 3 y) 



Z&> and m 3 = ^R e S ( 



ICO . 



00 = 



0o = 



sinh(m 3 ) 
sm(m 2 y) 



sin(m 2 ) 



"o = 



■ (e m ' (l + hm l ) - 1)+ - r - (hm 2 cotm 2 - 1) 



m 2 + m 2 



+ ■ 



. m 2 - m 2 



(/im 3 cothm 3 - 1) 



2 sinh mj - 2/imj coshmj 



2 (g'"' (l -hm l )-l)+ 2 ° r - (hm 2 cotm 2 - 1) 



m, 



+ ■ 



m t + m 2 



2 2 

y m l -m 3 



(hm 3 cothm 3 -l) 



-m,y 



2sinhm 1 - 2hm l coshm l 



X G r sin(m 2 _y) G m sinh(m 3 _y) 



+ — r + 



m 2 (m 2 + m 2 )sin(m 2 ) (m 2 - m 3 )sinh(m 3 ) 

G, 



m 2 + m\ 



[hm 1 cotm 2 -l) 



m 2 - m 3 2 



(/jffjj cothm 3 -l) 



sinhmj}' /i^coshm^y-^+sinhm^y-l) 



sinhwij -hm l coshm, 



( i A 



sinh m l - hn\ coshmj 



| X G r sm(m 2 y) | G ffl sinh(m 3 y) 
m 2 (mf + mjjsin^) (mf -m 3 )sinh(m 3 ) 



Again 



= 



sinh(m 3 ) 
m(m 2 y) 



sin 



( m 2 ) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



The dimensionless shearing stress at the upper wall of the channel is 
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A 

ml j 
G 



ml + ml 



{hm 2 cotm 2 -l) 



— (hm 3 cothm 3 -l) 



m l coshm. 



sinhmj -hm ] coshm, J (sinhm, -hm l coshm,)m, 



G r m 2 cot(m 2 ) G m m 3 coth(m 3 ) 

+ i 2 ; rr + r~2 rr - 

\m l +m 2 J \m l -m 3 J 



The dimensionless shearing stress at lower wall of the channel is 



m, 2 + ml 



(hm 2 cotm 2 - 1) 



-{hm 3 cothm 3 - 1) 



m. 



sinh m l - hm l coshm, 



N /l(coshmj - hm l sinhmj 
(sinh OTj - hm l coshm, )m l 



LI 



G r m 2 



G m m 3 



m, + ml 



)sin(m 2 ) (mf - m 3 )sinh(m 3 ) 
Rate of heat transfer across the channels at the upper wall is 



Nu x = 



m 2 cosm 2 iat 
sinm 2 



Rate of heat transfer across the channels at the lower wall is 



Nu 2 = 



dy 



m 



2 e U* 



J y=0 



sinm 2 



The rate of mass transfer at the upper wall is 

s\ = 



The rate of mass transfer at lower wall is 



m 3 coshm 3 to 



sinhm, 



Sh 2 = 



sinhm. 



e 



DISCUSSIONS AND GRAPHS 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



In this paper we have studied the effect of slip condition on unsteady MHD oscillatory flow in a channel filled 

with porous medium with heat radiation and mass transfer .The effect of the parameters Gr,h, Gm, H, N,w, D a , P r and S c 

on flow characteristics have been studied and shown by means of graphs and tables. To obtain the graphs and tables the 
modulii of velocity shearing stress taken w.r.t distance(y). Effect of different parameterson shearing stress at both the 
plates are shown on the same graphs. 

Figure-(l) illustrates the effect of the parameters Gr,Gm and H on velocity at any point of the fluid, when 
A.=l,Sc=0.22,S=l,h=l,w=l,t=l,Re=0.71,Pe=l and N=2. It is noticed that the velocity increases with the increase of Grash 
of number (Gr) and modified Grash of number(Gm) from the lower wall to the upper wall. Further the velocity at any point 
of the fluid decreases as the increase of Hartmann number (H) from the lower wall to the upper wall. 
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Figure-(2) illustrates the effect of the parameters X,Re and Pe on velocity at any point of the fluid, when 
H=l,Sc=0.22,S=l,h=l,w=l,t=l,Gr=l,Gm=l and N=2.It is noticed that the velocity decreases from the lower wall to the 
upper wall with the increase of both Renold number (Re) and peclet number (Pe). Further the velocity at any point of the 
fluid increases as the increase of the parameter (X) from the lower wall to the upper wall. 

Figure-(3) illustrates the effect of the parameters N,S and h on velocity at any point of the fluid, when 
w=l,Re=0.71,Gr=l,t=l,Gm=l,Pe=l,H=l, X =1 and Sc=0.22.It is noticed that the velocity at any point of the fluid increases 
with the increasing value of radiation parameter (N) from lower wall to upper wall. Further the velocity at any point of the 
fluid decreases with the increase of S and rarefaction parameter (h) from lower wall to upper wall. 

Figure-(4) illustrates the effect of the parameters w,t and Sc on velocity at any point of the fluid, when 
H=l,Re=0.71,S=l,h=l,w=l, X =l,Gi=l,Pe=l .Gm=l and N=2.It is noticed that the velocity decreases with the increase of 
frequency of oscillation(w) and Schmidt number(Sc) from the lower wall to upper wall. Further the velocity at any point of 
the fluid remain same with the increase of time (t) from the lower wall to the upper wall. 

Figure-(5) illustrates the effect of the parameters Gr, Gm and H on shearing stress at both upper and lower wall, 
when t=l, Re=0.71, Pe=l, Sc=0.22 h=l,w=l , S=1,X =1 and N=2.It is noticed that the shearing stress increases with the 
increase of Grashof number(Gr) and modified Grashof number(Gm), also decreases with the increase of Hartmann number 
(H) at both walls .The shearing stress at upper wall is less than lower wall. 

Figure-(6) illustrates the effect of the parameters Re, X and Pe on shearing stress at both upper and lower wall, 
when t=l, Gm=l, Gr=l,Sc=0.22 h=l,w=l, S=1H =1 and N=2.It is noticed that the shearing stress increases with the 
increase of parameter (X), also decreases with the increase of Renold number (Re) and peclet number (Pe) at both walls. 
The shearing stress at upper wall is less than lower wall. 

Figure-(7), illustrates the effect of the parameters S, h and N on shearing stress at both upper and lower wall, 
when t=l, Gm=l, Gr=l, Re=0.71 Sc=0.22, Pe=l, H =1 and w=l.It is noticed that the shearing stress decreases with the 
increase with of increasing S and rarefaction parameter (h), where as increases with the increase of radiation parameter (N) 
at both the walls. The shearing stress at upper wall is less than lower wall. 

Figure-(8) illustrates the effect of the parameters Sc and w on shearing stress at both upper and lower wall, when 
t=l, Gm=l, Gr=l, Re=0.71 h=l, Pe=l, S=l, H =1 and N=2. It is noticed that the shearing stress decreases with the increase 
both Schmidt number(Sc) and frequency of oscillation(w). The shearing stress at upper wall is less than lower wall. 




Figure 1: Effect of Gr, Gm and H on Velocity Profile when 
/-l, Sc=0.22, S=l, h=l, w=l, t=l, Re=0.71, Pe=l and N=2 
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Figure 2: Effect of Re, Pe and Ml) on Velocity Profile when 
H=l, Sc=0.22, S=l, h=l, w=l, t=l, Gr=l, Gm=l and N=2 
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Figure 3: Effect of N, h and S on Velocity Profile when 
w=l, Re=0.71,Gr=l, t=l, Gm=l, Pe=l, H=l, X =1 and Sc=0.22 
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Figure 4: Effect of Sc, w and t on Velocity Profile when 
H=l, Re=0.71, S=l, h=l, w=l, X =1, Gr=l, Pe=l. Gm=l and N=2 
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Figure 5: Effect of Gr, Gm and H on Shearing Stress when 
t=l, Re=0.71, Pe=l, Sc=0.22 h=l, w=l , S=l, X =1 and N=2 
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Figure 6: Effect of X (l),Re and Pe on Shearing Stress when 
t=l, Gm=l, Gr=l, Sc=0.22 h=l, w=l, S=1H =1 and N=2 
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Figure 7: Effect of N, S and h on Shearing Stress when 
t=l, Gm=l, Gr=l, Re=0.71 Sc=0.22, Pe=l, H =1 and w=l 
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Figure 8: Effect of Sc and w on Shearing Stress when 
t=l, Gm=l, Gr=l, Re=0.71 Sc=0.22, Pe=l, H =1 and w=l 



CONCLUSIONS 

This paper deals with the study of effect of slip condition on unsteady flow of an electrically conducting fluid 
trough a channel filled with saturated porous medium in the presence of transverse magnetic field, radiative heat transfer 
and mass transfer. Results are presented graphically to illustrate the variation of velocity, shearing stress, with various 
parameters. In this study, the following conclusions are set out: 

• The velocity increases for the increasing of Gr, Gm,N and X, also decreases for the increasing of Pe & Re ,Sc ,H,S 
and h from lower wall to upper wall. 

• The shearing stress of both plates increases with the increase of Gr, Gm,N and X, also decreases for the increasing 
of Pe & Re ,Sc ,H ,S and h. The shearing stress at upper wall is less than lower wall. 
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